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We discuss the question of gauge choice when analysing relativistic density perturbations at 
second order. We compare Newtonian and General Relativistic approaches. Some misconceptions 
in the recent literature are addressed. We show that the comoving-synchronous gauge is the unique 
gauge in General Relativity that corresponds to the Lagrangian frame and is entirely appropriate to 
describe the matter overdensity at second order. The comoving-synchronous gauge is the simplest 
gauge in which to describe Lagrangian bias at second order. 


I. INTRODUCTION 

Future galaxy surveys will probe scales comparable to the Hubble horizon, and it is therefore important to correctly 
incorporate relativistic effects in galaxy number counts. There are two types of relativistic corrections to Newtonian 
approximations. 

Firstly, relativistic corrections arise from nonlinear constraint equations in General Relativity (GR) [TH5] (see also 
®IZ! for related results). In a Newtonian analysis, the Poisson equation is a linear relation between the gravitational 
potential and the matter overdensity. In GR, this is replaced by a nonlinear relation between the metric perturbation 
and the matter overdensity. Due to this nonlinear relation, the density field becomes non-Gaussian even if the 
primordial curvature perturbation is Gaussian. The non-Gaussianity in the comoving-synchronous density field can 
be parametried by the effective local parameter 00 

/is. R =-|. a) 

The second type of GR correction arises because we observe galaxies on the past lightcone and not a constant 
time hypersurface. On Hubble scales, the number overdensity is gauge-dependent, and none of the standard gauge- 
invariant choices corresponds to the actual observed overdensity, which is automatically gauge-invariant. This observed 
overdensity includes all redshift, lensing and volume distortions, some of which make a growing contribution on large 
scales [8HT0] . The first-order GR lightcone effects contaminate the signal from primordial non-Gaussianity and so 
these effects must be included in order to extract the primordial signal mm- 

Recently, several papers have extended the calculation of the observed galaxy number overdensity on cosmological 
scales up to second order in relativistic perturbation theory [171421] . These results could be important to make accurate 
estimates of cosmological parameters, including primordial non-Gaussianity, from large-scale structure. 

Since the observed galaxy overdensity does not depend on gauge, we can calculate it using any gauge. However, 
a subtlety arises when we define the galaxy bias. Local bias should be defined in the rest-frame of cold dark matter 
(CDM), which is assumed to coincide with the rest frame of galaxies on large scales. Local bias can be computed 
using the peak-background split approach, where halo collapse occurs when small-scale peaks in the density exceed 
a critical value [ 22 ] . Long-wavelength modes modulate this critical value. The spherical collapse model has an exact 
GR interpretation and the criterion for collapse of a local overdensity is when the linearly evolved density in the 
comoving-synchronous gauge reaches a critical value |23H26| . Thus it is natural to define bias in terms of the density 
in the comoving-synchronous gauge mm- At first order, the galaxy overdensity in comoving-synchronous gauge is 
given on large scales and for Gaussian primordial fluctuations by 

S g c(z,x) = b(z)5 c {z, x), 

where Sc is the matter overdensity. This approach was extended to second order in PH- 
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Recently, the choice of the comoving-synchronous gauge at second order has been criticised f2TU29] . with claims 
that it leads to violation of mass conservation and is inappropriate for defining galaxy bias. Instead, [27] advocated 
that the correct gauge to choose is the total matter gauge [25H5U] . Here we show that, on the contrary, 

• there is no violation of mass conservation in comoving-synchronous gauge, 

• the comoving-synchronous gauge is the appropriate gauge for defining the local Lagrangian bias at second order. 
We also clarify the similarities and differences between the Newtonian and GR approaches. 


II. NEWTONIAN FRAMES VERSUS GR GAUGES 

In this section, we first define the Eulerian and Lagrangian frames in Newtonian theory and the corresponding 
comoving-synchronous gauge in GR. Then we consider a transformation to a gauge in GR that corresponds to an 
Eulerian frame. The issue of residual gauge freedom in the comoving-synchronous gauge is also discussed. 

Throughout we assume that the matter is in the form of irrotational dust, representing CDM. This is a reasonable 
assumption on cosmological scales. We use conformal time r in a Friedmann-Lemaitre-Robertson-Walker background 
spacetime, with scale factor a(r) and H (r) = a'(r)/a(r). 


A. LAGRANGIAN FRAME AND COMOVING-SYNCHRONOUS GAUGE 


Newtonian theory 


In Newtonian theory, the equations that govern the dynamics of irrotational dust in an Eulerian frame, (r, x l ), on 
an expanding background are 


7 C 

—+ (1 + Sn)V 1 vni = 0 (continuity), 
dr 

+ Hv l N + V V = 0 (Euler), 

dr 

V 2 </> = 47t Gpa 2 SN (Poisson), 

where d/dr = d/dr + v l N V i is the convective derivative and the velocity is derived from a potential: 


In terms of the deformation tensor, 


■ l3k diV NJ =0 => v l N = d l v N . 


= d l v J N = + a%, 


( 3 ) 
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where d is the perturbation to the expansion and er® is the shear, the continuity and Euler equations can be written 


as 


d$N 

dr 


+ (1 + <5/v)$jv — 0 , 


d^N 

dr 


+ 'H'disr + + AnGp6]s[ — 0 , 


( 8 ) 
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where we used the Poisson equation. 

In the Lagrangian picture, the dynamics is described by comoving coordinates q l , specified at an initial time, so 
that each fluid element maintains the same coordinates during evolution. Eulerian coordinates x l are related to the 
Lagrangian q l by 


r i (q,r)=g i + ^(q,r), 


( 10 ) 


where the displacement field 'L® vanishes at the initial time, and determines the velocity in the Eulerian reference 
frame: 


Vn = v 1 e = 


d 4®® 

Ih 


( 11 ) 
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The matter density is a scalar quantity that can be written as a function of either Eulerian coordinates or Lagrangian 
coordinates. These functions will have different functional dependence on the corresponding coordinates, but the 
density is the same physical quantity at a given point. Hence we can write the density contrast at any fixed point as 

5 n = = <Mx(q,T),r). (12) 


GR Cosmology 


It is possible to define a Lagrangian frame uniquely for irrotational dust in GR 1 . This is given by the comoving- 
synchronous gauge, with comoving coordinates and constant time hypersurfaces r = const orthogonal to the fluid 
worldlines. The line element is 


ds 2 = a(r) 2 [ - dr 2 + 7 ^ (r, q k )dq l dq J ] = o(t) 2 [ - dr 2 + exp (2([(t, q k )) %j{r,q l )dq l dq J ], (13) 

where £ is the comoving curvature perturbation and ( 7 ! = 1 . 

The fluid four-velocity is given by u M = — aS ^, which coincides with the unit normal to r = const. The deformation 
tensor is d£ = aV^zt^ — n(5£ + zz^zz,,), which in the comoving-synchronous gauge is purely spatial and is given by the 
extrinsic curvature [31j 




(14) 


where a prime denotes d/dr. The Einstein equations in this gauge give 

z ? 2 — d x h+ 4 Hd + R ^ = 167 rGa 2 p< 5 , (energy constraint), 
D J dij = did, (momentum constraint), 

while the evolution of the deformation tensor is given by 

d )' + 2Hd) + dd) + 1 (d kl d kl - d 2 ) 8) + Rf ]i - ^R^S) = 0, 
?(3) 


(15) 

(16) 


(17) 




where R\j is the Ricci tensor of 7 ^ and ID,; is the covariant derivative defined by 7 ^. By combining the trace of ( |17[ ) 
and the energy constraint (15), we obtain the Raychaudhuri equation 


d' + Rd + dt’du + 47rGa V = 0 . 


(18) 


The continuity equation is 


6' + (1 + S)d = 0 . 


(19) 


The two equations (18) and (19) are formally equivalent to the Newtonian equations in Lagrangian form, if we replace 
the comoving-synchronous time derivative in GR, d/dr, with the convective derivative in Newtonian theory, d/dr. 


B. EULERIAN FRAME AND TOTAL MATTER GAUGE 


In GR, any frame can legitimately be called Eulerian, while there is a unique Lagrangian frame for irrotational 
dust, i.e. the comoving-synchronous gauge as discussed above (and see [22). However, it is possible to define a gauge 
where the density contrast obeys the same evolution equations as the Newtonian theory in the Eulerian frame up 
to the second order (3[ {28, 33] ■ We follow [34] [353 i n calling this gauge the total matter gauge. It is also called the 
velocity orthogonal isotropic gauge [ 21 ] and the temporal comoving spatial G-gauge m- 

In Newtonian theory the transformation between the Eulerian and Lagrangian frames is given by a spatial change 
of coordinates, (10). We can apply the same spatial transformation to the comoving-synchronous gauge in GR, 


and this induces a velocity perturbation in the total matter gauge through (11), with d&T = d^i/dr. Since (10) 


1 It is possible to choose an irrotational growing mode for dust at first order, but at second order dust does not remain irrotational. 
Nonetheless we can still impose the comoving-orthogonal gauge conditions at second order with respect to the irrotational part of the 
velocity field. 
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represents a purely spatial gauge transformation in GR, the two gauges share the same time-slicing, i.e. the constant 
time hypersurfaces r = const are the same in the two gauges. 

There are in general two scalar degrees of freedom in the metric at each order H3 [351 13ZI- In the comoving- 
synchronous gauge these are both part of the spatial metric, (13), which (ignoring vector and tensor modes) can be 
written as 


7 Cij = (1 + 2 i>c)8ij + 2 Ec,i 


( 20 ) 


The gauge transformation (10) induces a shift (goi) in the metric, and can be used to make the spatial metric 

( 21 ) 


conformally flat, so that the total matter gauge is completely fixed by 

4>t = 4>c, B t = v T , E t = E c — J dr'u T (x, r') = 0, 

and the line element in the total matter gauge at first order becomes 2 

ds 2 = a 2 [ — dr 2 — 2 diB T dx l dr + (1 + 2 ip T )S i jdx l dx : ’] . 

The density contrast and the trace of the deformation tensor (representing the perturbation in the expansion) do not 
change under a purely spatial gauge transformation at first order: 


( 22 ) 


S T = Sc, At = tic- 


(23) 


Note however that in the comoving-synchronous gauge we identify i9c = \7 2 E' C , while the total matter gauge is fixed 
by setting Et = 0 and instead we find dr = V 2 ut- 

The total matter gauge fixes the gauge completely but there is a residual gauge freedom in the spatial gauge 
transformation to the comoving-synchronous gauge, 


/ T pT 

dfvr{x, f) = / dfvx{x, t) + e(x). 


(24) 


This can be fixed by choosing the spatial coordinates q l in the comoving-synchronous gauge to coincide with the 
Eulerian spatial coordinates x l (i.e., the spatial coordinates in the total matter gauge) on an initial hypersurface r ln = 
const., so that 


e(x) = 0. 


(25) 


In other words, by an appropriate choice of initial conditions, this gauge mode in the comoving-synchronous gauge can 
be removed 3 . This is equivalent to the usual Newtonian assumption that the Eulerian and Lagrangian coordinates 
coincide at an initial time, and hence the displacement field is zero initially, \IT(q, r; n ) = 0. In particular, to study 
the growing mode of primordial density perturbations (the regular mode at r- ln —> 0), we will require e(x) = 0 in the 
limit r; n —> 0. 

Just as in Newtonian gravity, the density perturbation is a 3-scalar so remains invariant at a fixed physical point 
on the constant-time hyper-surface under the spatial coordinate transformation (10) 


At second order this gives 


Mq,T) = Mx(q,r),r) • 

J T (x, r) = Sc(x,r) - ^(x,r)diSc(x,T) 

= S c (x,t)+ [d. l S c (x,T)]V~ 2 d z S c (T,x), 


where in the last line we have used (111 and the continuity equation (19) at first order, 

$' + 0 = 0 , 


(26) 


(27) 


(28) 


to eliminate the displacement field in favour of the density contrast. It is clear that the resulting transformation 
(271 is highly nonlocal at nonlinear orders. 


2 In this work, we define the total matter gauge at nonlinear order as the gauge that shares the same constant time hypersurfaces as the 
comoving-synchronous gauge (i.e., normal to the comoving worldlines) but with spatial coordinates such that Et = 0 at all perturbation 
orders. 

3 To quote Bertschinger ESI: “The presence of these extraneous solutions (called gauge modes) has created a great deal of confusion in 
the past, which might have been avoided had more cosmologists read the paper of Lifshitz ESI”. 
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C. NEWTONIAN AND GENERAL RELATIVISTIC CONSTRAINTS 

In Newtonian theory, the evolution equation for is obtained by taking the spatial derivative of the Euler 
equation, which contains the second spatial derivative of the Newtonian potential, the tidal tensor. There is no 
evolution equation for the tidal tensor in Newtonian theory and it is thus impossible to obtain a complete system of 
evolution equations. This is due to the underlying elliptic nature of the system, i.e. it embodies action-at-a-distance. 
The system is closed by a constraint, i.e. the Poisson equation. Various approximations close the evolution system 
using some suitable assumptions. For instance, in the Zel’dovich approximation m the evolution equation for the 
tidal tensor becomes redundant, and the continuity and Raychaudhuri equations, together with the evolution equation 
for the shear, form a closed system of differential equations (see e.g. |H]). (For a related discussion on obtaining the 
Newtonian equation from the GR one, see [43] .) 

In GR, a covariant approach can be used, via a 1+3 split based on the fluid four-velocity HUES)- This leads 
to a system of evolution and constraint equations for the kinematical variables and the electric and magnetic parts 
of the Weyl tensor, E^ and H^. The whole set of equations is hyperbolic, i.e. it embodies causality. In GR the 
evolution equations for and close the system of evolution equations. 

Another important difference lies in the constraint equations. In the Newtonian case, the only constraint is the 
Poisson equation, which provides a linear relation between overdensity and potential. In GR, there are two constraints, 
the energy constraint (151 and the momentum constraint (16). At first order in perturbations, the two constraints in 


GR combine to give the relativistic Poisson equation, 


V 2 4- = 4tt Ga 2 p5 c , 


(29) 


where 4* is the metric perturbation in longitudinal gauge (see e.g. [35 ] 145 ] for the gauge-invariant version). 

However, a difference arises at higher orders W- At higher orders, the density field obeys second order differential 
equations that are sourced by the product of the lower order quantities. The time evolution of the homogeneous 
solution is given by the linear growth function D + (t) and its momentum dependence is determined by the constraint 
equation. The particular solution on the other hand describes the nonlinear evolution of the initial density contrast. 
The nonlinear density contrast in comoving-synchronous gauge may be written as 


Mq, T ) = C(q)D + (r) + Pn(<l)D n+ (q, r), 


(30) 


n =2 


where C(q)D + (r) is the growing mode solution of the homogeneous part of the evolution equation for the density 
contrast at all orders and the sum represents the nonlinear evolution of the initial density contrast at any given point 
q in comoving coordinates. 

On large scales, the growing mode amplitude C(q) is related to the curvature perturbation £ by the energy constraint 
equation (15) @|: 


C(q) oc exp(-2C) - 4V 2 C - 2(V£) 


(31) 


Even in the absence of primordial non-Gaussianity, the higher-order homogeneous solution is not zero. At second 
order, this leads to the effective non-Gaussianity of local type in 0 nHg. By contrast, in the Newtonian 
approximation, the Poisson equation is a linear relation between the density contrast and Newton potential at all 
orders, and therefore the higher-order homogeneous solution is zero in the absence of primordial non-Gaussianity. 


III. IS THE COMOVING-SYNCHRONOUS GAUGE A BAD GAUGE? 

The comoving-synchronous gauge has been described as inappropriate to use for analysing the growth of large-scale 
structure at nonlinear order fTTU'IU] . The specific claims include: 

• comoving-synchronous gauge “violates mass conservation” ED; 

• the density contrast in comoving-synchronous gauge cannot be used for galaxy bias mm- 
In this section, we show that these claims are not correct. 
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A. MASS CONSERVATION 


Newtonian theory 

In Newtonian theory, the spatial average of the overdensity in Eulerian coordinates is mm 

(Se(x,t)) x = ^ d 3 xS E (x, t) , 

where the average is made over a volume that goes to infinity. The Lagrangian spatial average is 

{S L {q,r)) q = ^ J d 3 q J - 1 S L (q,T ), 

where we must take account of the Jacobian for Lagrangian coordinates 

Bx i 

J = 


dq'i 


( 32 ) 


(33) 


(34) 


Taking into account that £e(x(q, r), r) = Sl{ q, t), we find the relation between Eulerian and Lagrangian averages: 

(S e (x,t)) x = {S L (q,T)) q . (35) 

GR cosmology 

In GR, we use T- and C-gauges to correspond to the Eulerian and Lagrangian frames, as discussed above. The 
Eulerian average is written using the covariant volume element: 


(6 t (x,t)) x = ^ d 3 x |7 t| 1/2 $t(x, t), 

where 7 t is the metric induced on V in the T-gauge, while in the Lagrangian C-gauge we have 

(Mq,-r)) g = ^ d 3 q \lc \ 1/2 fc(R r). 

Note that at first order, 

l7d' /2 


175 


11/2 


= 1 + V Eq ■ 


and hence using (241 and we can identify this with the Jacobian ( |34[ ) 

J = \ + d i * i . 


(36) 


(37) 


(38) 


(39) 


Actually, J = (|7c|/|7t|) 1//2 is a general result in tensor calculus for the coordinate transformation of the determinant 
of a metric 7 ^. In the transformation to the comoving-synchronous gauge, the volume element is invariant 


d 3 x | 7 t | 1/2 = d 3 q | 7 c | 1/2 , 

since the change of coordinates is purely spatial. Thus we have, as in Newtonian theory, 

(S t (x,t)) x = (S c (q,r)) q . 


(40) 


(41) 


The spatial average of the density contrast in the total matter gauge must equal the spatial average in the comoving- 
synchronous gauge since we are evaluating the same physical quantity on the same physical hypersurface, just using 
different spatial coordinates. 

The solution up to second order in the total matter gauge (restricting for simplicity to the case of planar symmetry 
or matter-dominated era) can be written in Fourier space as 

S T (k,r) = D+(r)S (1) (k) + Dl(r) J d 3 kF 2 T (k, k - k) J ( 1 ) (k) jW(k-k). 


(42) 
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Here we use the Fourier transform with respect to Eulerian coordinates 

d 3 x 


m = ] 

The Eulerian kernel vanishes in the large-scale limit 


„ — zk-: 


(2tt) 3 


: s(x) ■ 


lim F 2 T (ki,k 2 -ki) = - 

k 2 —>0 I 


3-5 


(ki ■ k- 

h 2 h 2 


™ k? 


2 

kV 


(43) 


(44) 


and hence the spatial average vanishes in the large-scale limit. 

In [23 it is argued that the kernel for the comoving-synchronous gauge density contrast does not vanish in the 
large-scale limit: 


lirn F 2 c (ki, k 2 - k x ) = 1 + - 
&2—>-0 ( 


(ki • k 2 ) 2 
h 2 lc 2 


-1 


k\ 

k\ 


(45) 


and thus the spatial average does not vanish at second order, violating mass conservation. 

This misunderstanding comes from taking the Fourier transform in Eulerian coordinates, corresponding to the 
spatial average (<5c(x, T )) x - Using the second-order gauge transformation |27| ), we note that 

(46) 


(Mx,r)) x = (Mx,t) - ’F l (x,T)d i £c(x,T)) x = (<5c(q.r)) 


Thus using (11) and (|28j), we find that the apparent difference in the spatial average of St and Sc in Eulerian 


coordinates is due to a dipole term, which can be written in Fourier space as 


['F i 3i<5c](k,T) = D 2 + {t) 


d 3 ki 

(27r) 3 


kl k 



fc(ki)<5c(k 2 ) , 


(47) 


where k 2 = k — k^. Clearly this does not vanish in the large-scale limit, k —> 0. Because of this dipole term, )27| 
concludes that the matter density fluctuation in the comoving-synchronous gauge violates mass conservation at second 
order. 


The problem with naively using (45) in (36) is that this is the Eulerian average relation, instead of the correct 


Lagrangian average (371. The result (47) shows clearly that the dipole term is not missing in comoving-synchronous 


gauge - we can recover the same kernel that we get in the total matter gauge (see, for example, the appendix of 0 )- 
This proves that, in both the total matter and the comoving-synchronous gauges, the kernel F 2 vanishes as k 2 in the 
large-scale limit. Therefore we recover the correct result in either gauge without any violation of mass conservation 
(which is of course respected in GR in any gauge). 


B. LAGRANGIAN GALAXY BIAS 

In Newtonian theory there are two common prescriptions for describing the distribution of galaxies (assumed to 
reside in collapsed dark matter halos) as biased tracers of the underlying matter density. An Eulerian bias model 
assumes that the local overdensity of galaxies traces the local matter overdensity at a given Eulerian point x, while 
in a Lagrangian bias model it is the initial overdensity in the Lagrangian coordinates that determines the galaxy 
distribution. 

As we have seen, the total matter gauge in GR corresponds to an Eulerian frame. On the other hand, Lagrangian bias 
models follow the collapse of initial overdensities in the matter field at a fixed point q in Lagrangian coordinates, i.e., 
comoving with the matter distribution (which moves through Eulerian space). Therefore the comoving-synchronous 
gauge provides the simplest GR gauge in which to describe Lagrangian bias. In the Press-Schecter theory [22], for 
example, halos of a given mass, M, collapse when the linearly growing local density contrast <5 lln (q) (smoothed on the 
corresponding mass scale) reaches a critical value. In particular, spherical collapse has an exact GR solution in the 
comoving-synchronous gauge [25] . 

The local galaxy number density in comoving-synchronous coordinates, n g c{ q), is thus a function of the initial 
density field, corresponding to a local Lagrangian bias model [461449] 

^pc(q) = &i'6 lin (q) + & 2 [<5 lin (q)] 2 + ... . 


(48) 




















Since the initial displacement between Lagrangian and Eulerian frames is zero, at sufficiently early times where we can 
also neglect nonlinear evolution, the density held is the same in both comoving-synchronous and total matter gauges, 
where x = q initially. However as the density contrast evolves, a local Lagrangian bias model in general predicts a 
nonlocal Eulerian bias mm- 

Note that the Lagrangian bias model describes the number density n g with respect to an initial (Lagrangian) 
coordinate volume and thus n g is not a 3-scalar in GR, but rather a density of weight one. Transforming the number 
density from Lagrangian to Eulerian coordinates, we must include the Jacobian: 

n gT (x) = n gC ( (49) 

and thus in terms of the galaxy overdensity we have 


1 + 6 gT (x) = [1 + <5t(x)][1 + <5 gC (q(x))]. 


(50) 


Transforming the Lagrangian galaxy density contrast 5 g q from Lagrangian q to Eulerian x using (271, then introduces 
nonlocal terms at second and higher order in the Eulerian galaxy density contrast <5 9 t(x, t), as shown for example in 


6 gT = b?6 T + b$[6 T ] 2 + b s s 2 + ..., (51) 

where s 2 describes second-order tidal terms. 

The local Lagrangian bias model has been shown to give a better fit to the distribution of collapsed halos in N-body 
simulations than a local Eulerian bias model (5T1 [52] . For complete generality, the Lagrangian bias model may also 
include nonlocal terms due to tidal effects [53[. 

The initial density contrast (in Lagrangian or Eulerian frames) is dominated by the linearly growing mode 4 : 

S hn (q, r) = C(q)D + (r). (52) 


In GR, unlike in Newtonian gravity, the growing mode amplitude C(q) is nonlinearly related to the primordial 
curvature perturbation £ by (31). Thus even if £ is Gaussian, the initial density field is non-Gaussian. There has been 
significant interest in the effect of non-Gaussianity of the initial density field in GR on galaxy bias m- It is beyond 
the scope of this paper to discuss in details this effect and we leave it to future work. 


IV. CONCLUSIONS 

Recently, several papers have extended the calculation of the observed galaxy number overdensity on cosmological 
scales up to second order in relativistic perturbation theory H2JSII- These results include all GR effects that arise 
from observing on the past lightcone and will be important for accurate cosmological parameter estimation, including 
non-Gaussianity. In order to pursue such applications, it is important to carefully choose the gauge in which to relate 
the fluctuations of galaxy number density to the underlying matter density fluctuation. 

For an irrotational dust flow in ACDM, assuming that the galaxy velocity is equal to the CDM velocity on large 
scales, El claimed that a suitable gauge to define the local bias up to second order is the comoving-synchronous 
gauge. This choice was criticised in EZj, with the claim that the comoving-synchronous gauge “does not properly 
represent the matter density fluctuation with the mean at the local proper time, violating mass conservation”, and 
the further claim that the correct gauge to choose is the total matter gauge. 

We have shown that the comoving-synchronous gauge: 

• is the gauge in GR that best matches the concept of Lagrangian frame in Newtonian theory, 

• does not violate any mass conservation at second order, and 

• is appropriate for defining local Lagrangian galaxy bias up to second order. 

As part of our argument, we have highlighted the differences between the Newtonian approximation and the correct 
General Relativistic analysis. 


4 Here the “linearly growing mode” indicates that the time evolution is given by the linear growth function D- |_(t). 
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